Spontaneous autophoretic motion of isotropic particles 
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Suspended colloidal particles interacting chemically with a solute are able to sell-propel by au- 
tophoretic motion when they are asymmetrically patterned (Janus colloids). Here we demonstrate 
that the chemical anisotropy is not a necessary condition to achieve locomotion. The non linear 
interplay between surface osmotic flows and solute advection can produce spontaneous, and sell- 
sustained motion ol isotropic particles. We solve, lor a spherical particle, the classical nonlinear 
autophoretic theoretical framework at arbitrary Peclet number. For a given set ol material pa- 
rameters, there exists a critical particle size, or equivalently a critical Peclet number, above which 
spontaneous autophoretic motion occurs. The flow induced by the particle lurther displays a hier- 
archy ol instabilities associated with quantized critical Peclet numbers. Using numerical solutions 
ol the lull (unsteady) diflusiophoretic problem we confirm our analytical predictions and show that, 
above the instability threshold, the isotropic particles reach a steady swimming state with bro- 
ken front-back symmetry in the concentration field and the hydrodynamic signature ol a "pusher" 
swimmer. This instability to propulsion could be relevant to the high-throughput production ol 
sell-propelled particles. 

PACS numbers: 82.70.Dd, 47.70.Fw, 47.20.-k, 87.85.gj 
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The locomotion of microorganisms has long been used 
as a motivation and a practical inspiration for the design 
of synthetic self-propelled particles. Typically, biological 
cells generate propulsion by deforming their slender ap- 
pendages, termed fiagella or cilia, in a non-time-reversible 
fashion However, and perhaps not surprisingly given 
the numerous microfabrication challenges, no genuinely 
self-propelled micro-swimmer has been manufactured in 
the lab so far. Instead, man-made biomimetic propellers 
are driven by external torques or forces. That actuation 
allows either to deform softpropellers whose deformed 
shape induce propulsion [2j-|4j, to continuously generate 
propulsion in chiral shapes [a, Q , or to exploit interac- 
tions with surfaces 

An alternative route for the production of artificial 
small-scale swimmers has proven to be much more suc- 
cessful. It consists in making miniaturized chemically 
powered "engines" with no movin g pa rts, typically made 
of reactive Janus beads or rods |10H12| . The reaction 
products released by these chemically-asymmetric parti- 
cles create concentration gradients which induce a net 
phoretic fluid motion near their surface leading to loco- 
motion. Theoretically, the interest in these so-called au- 
tophoretic swimmers was triggered by a theoretical model 
which accounted for such a novel propulsion mechanism 
in a simple and generic fashion [13J. This model was 
then further elaborated to include the nonlinear inter- 
play between the colloid motion and the advection of the 
reactants 14, to deal with the rotational Brownian 



motion of the swimmers [ly], and to detail the micro- 
scopic coupling between the concentration gradients and 
the fluid flows at small scales 
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In order to self-propel, autophoretic swimmers are 
chemically patterned, and it is the asymmetries in the 
chemical reactions on the swimmer surface which are 
responsible for locomotion in the first place. This re- 
quirement would make it thus difficult to achieve high- 
throughput production. An ingenious solution to such an 
engineering issue was recently offered with the produc- 
tion of isotropic self-propelled Marangoni droplets 18 1. 



In a mechanism akin to the one responsible for the spon- 
taneous motion of reactive droplets surfing on fluid inter- 
faces [lsl, H(| , a net flow is induced around the droplets by 
interfacial stresses induced by self-generated gradients of 
reactive surfactants [H , [21 1 . Unlike autophoretic swim- 



mers, which swim as a result of built-in asymmetries in 
their design, the concentration of surfactant molecules at 
the surface of reactive droplets is spontaneously broken, 
leading to propulsion. 

In this letter we demonstrate that autophorcsis can 
produce spontaneous motion of isotropic solid particles 
even in the absence of interfacial (Marangoni) stresses. 
We solve analytically the classical nonlinear autophoretic 
theoretical framework at arbitrary Peclet number for per- 
turbation around the isotropic, purely diffusive, state. 
We show that, for a given set of material properties, there 
exists a critical particle size above which spontaneous au- 
tophoretic motion occurs. In addition, the flow induced 
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by the reactive particles displays a hierarchy of instabil- 
ities associated with quantized critical Peclet numbers. 
Using numerical solutions of the full unsteady diffusio- 
phoretic problem we confirm our analytical predictions 
and show that, above the instability threshold, isotropic 
particles reach, after a transient, a steady swimming 
state with broken front-back symmetry in the concentra- 
tion field and the hydrodynamic signature of a "pusher" 
swimmer. 

We consider a spherical particle of radius a located in 
a Newtonian fluid of shear viscosity rj. We focus on the 
limit of steady Stokes flow so that the inertia of both 
the fluid and the particle is negligible. The sphere is 
chemically active and cither emits or captures a single 
type of solute particles with an isotropic surface emission 
rate denoted A ljj l22[. Within the standard minimal 
framework, we ignore in what follows the specifics of the 
chemical reaction occurring at its surface, and all other 
reactants and reaction products, which are assumed to 
interact only weakly with the particle. The solute in- 
teracts with the spherical particle through a short-range 
potential. The competition between the interaction po- 
tential and molecular diffusion sets the thickness of the 
so-called diffuse layer, A. As in classical approaches, we 
consider the thin-diffuse layer approximation, A <C a [23f . 

Tangential gradients in solute concentration induce a 
net slip (phoretic) velocity outside the diffuse layer, to- 
gether with a tangential (Marangoni) stress discontinu- 
ity at the particle surface [II, [23( ■ For a solid particle, 
Marangoni stresses vanish, and the phoretic slip velocity 
is set by the local chemical-potential gradient parallel to 
the surface j24[. For small variations of the concentra- 
tion, c(r,t), its azimuthal component is given by 



, . M dc 

ue{r = a) = -w 



(1) 



The mobility, Ai ~ ifc^TA 2 /rj, can either be posi- 
tive or negative depending of the form of the solute- 
surface interactions (kn is the Boltzmann constant, and 
T the temperature) [23(. Assuming that the solute has 
a molecular diffusivity T>, the typical autophoretic ve- 
locity which sets the swimming speed of a Janus colloid 
is V = | AM [22j]. Non-dimensionalizing velocities, 
lengths, and concentrations in the transport equations 
by V, a, and a\A\/T> respectively, the coupled fluid flow- 
solute advection-diffusion problem takes the form 



V 2 u = Vp, 



V • u = 0, 

V 2 c, 



(2) 
(3) 



where we have defined the (signed) Peclet number as 
Pe = AM-a/V 2 . We henceforth consider axisymmetric 
solutions only, so that in spherical polar coordinates the 
solute concentration is written as c(r, /i, t) with fi = cos 8, 
and only rectilinear motion along the axis of symmetry 



e 2 is considered. In dimensionlcss units, the surface ac- 
tivity and the particle mobility become unitary numbers, 
A = A/\A\ = ±1, M = M/\M\ = ±1, and the boundary 
conditions on the particle surface become 



u(?' — > oo) = — Ue z , c(r 
dc 

u r (r = 1) = 0, 



■-(r = l) = -A, 



-dc 



ug{r = 1) = -My/l-[j? — (r = 1). 



(4) 
(5) 
(0) 
(7) 



In Eq. Q , U(t) is the dimensionlcss swimming velocity, 
obtained by enforcing that the particle is force-free. Us- 
ing the reciprocal theorem, U (t ) is found to be the surface 
average of the slip velocity [H, [IE] and, using Eq. ([T]), is 
written in terms of the first moment of c(l, /i, t) only as 



U(t) = —M / fj,c(l,fx,t)dfj,. 



(8) 



The dimensionless problem is then fully characterized by 
the signs of both A and M, the far-field concentration, 
Coo, and the value of |Pe|. 

In the case of uniform surface activity (constant value 
of A), a trivial solution exists at all Pe numbers, namely 
the solute concentration is isotropic, c = A/r + Coo, lead- 
ing to no net flow (u = 0) and zero swimming velocity 
(U = 0). However, both the solute transport equation, 
Eq. ©, and the boundary conditions on the particle, 
Eq. J7|), couple the swimming problem with the solute 
dynamics. A small fluctuation of the particle velocity 
would result in a polar perturbation of the flow field. 
Due to the nonlinear convective coupling, u • Vc, this ve- 
locity fluctuation would lead to a polarization of the con- 
centration field around the finite-size particle. In turn, 
the first moment of the surface concentration, c(l,fj,,t), 
would become finite and, depending on its sign in relation 
to the initial perturbation, increase or decrease the par- 
ticle velocity through Eq. (|8]). If the velocity decreased, 
the initial perturbation would be stabilized, and no net 
motion could occur as a result of an infinitesimal fluc- 
tuation. However, if the velocity increased, the broken 
symmetry in the solute concentration would be amplified, 
and spontaneous motion would occur. 

To quantify the conditions for spontaneous motion, we 
analytically investigate the stability of the isotropic state. 
Defining c = c + c', u = u', U = U', and subsequently 
dropping the primes to denote perturbations, the Stokes 
flow problem around the sphere can be solved analyti- 
cally using the decomposition of the slip velocity along 
the so-called squirming modes [2(| 27 [. We then project 



both the advection-diffusion equation and the definition 
of the slip velocity along these modes, and take the first 
moment of the resulting equation with respect to the an- 
gular dependence, fi. This leads to an inhomogeneous 
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eigenvalue problem for the first moment of the concen- 
tration profile, ci(r, t) = | /uc(r, ^, t) d/^, as 



|Pe 



gci 1 
Si r 2 



, <9ci 
9r 



2 Cl 



J7Pc / 1 



U = --Mc x (l), 



Ci (l) = 0. 



(9) 
(10) 



Note from Eq. §E§ that the first moment, ci, evaluated 
at r = 1 is proportional to the instantaneous swimming 
speed of the particle, (Eq.Qj)]). Looking for eigenmodes of 
the form ci(r, t) = e cr *c(r), it can be shown that all eigen- 
values a of Eqs. (|9])- ([T0| are real, and that any a < is 
a solution. We focus here exclusively on potential unsta- 
ble modes (a > 0) and define ft = 'tr|Pc| > 0. Intro- 
ducing the rcscaled radial variable x = /3r, the function 
C(x) = c(r/f3) satisfies 



d 

d.r 



,dC 
d.r 



(2 + x 2 ) C 



2Pc 
~3~ 



1 



(11) 



The general solution of Eq. (fTTj) satisfying the far-field 
condition C(x — > oo) = is 



be~ x (l + x) 



A{x) B{x) 
8x 2 8x 



1 

4x 3 



(12) 
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FIG. 1: (Color online) Growth rate of the unstable swimming 
mode (theoretical prediction, solid line) as a function of the 
Peclet number, Pe, for emitting particles (A = 1) and pos- 
itive mobility (M — 1). Spontaneous symmetry-breaking of 
the concentration field and swimming occur at Pe = 4. The 
crosses represent the growth rate of the swimming mode n — 1 
as obtained from numerical simulations of the full unsteady 
coupled problem. The growth rate of the unstable modes of 
azimuthal order n — 2, 3 and 4 are also shown (theoretical 
prediction, dashed, dash-dotted and thin-solid, respectively). 
Note that the hierarchy in the growth rates is reversed at high 
Pe. 



where b is an integration constant to be determined, 
and A(x) = sinh(.T)Chi(x) — cosh(a;)Shi(2;), and B(x) = 
sinh(x)Shi(a;) — cosh(x)Chi(.T), with Chi(a;) and Shi(.T) 
being the hyperbolic cosine and sine integral functions 
respectively [28[ . Applying the two boundary conditions 
on the sphere, C"(/3) = 0, and C(/3) = 1, and using the 
definitions of A[x) and B(x) yields a final implicit ex- 
pression for the growth rate, a = /3 2 /|Pe|, as a function 
of the signed Peclct number 



Pe 



12/3 2 + 24/3 + 24 



/3 J 



(13) 



dt + 6 + f3 2 - /3 3 - 2(3 



For positive values of 0, the right-hand side of Eq. (TIB")) 
is strictly greater than 4. For any value of Pe above this 
critical value, a single positive value of (5 exists such that 
Eq. (|13p is satisfied. Consequently, the fluctuations of the 
first moment, ci, and the particle swimming speed, U , are 
exponentially amplified if the condition Pe > Pei = 4 is 
satisfied. Given that Pe is a signed quantity, the insta- 
bility condition requires that MA = 1. Particles with 
positive (resp. negative) mobility M. are unstable only if 
they have a positive (resp. negative) flux A correspond- 
ing to the case of emitting (resp. absorbing) the solute 
on the particle surface. The growth rate of the unstable 
swimming mode is shown as a function of Pe in Fig. [TJ 



(solid line). From a practical point, the existence of crit- 
ical Peclet number implies that there exists a critical 
particle radius, a\ = 425/V, above which an isotropic 
reactive particle at rest would undergo spontaneous mo- 
tion. For instance, considering the prototypal example 
of a Pt-coated bead in a H2O/H2O2 solution, we have 



V « 20 /ims-' [nll2ll29| and taking V w 1(T 9 m 2 /s [30| 
we obtain a critical radius a\ w 100 /iin, which is experi- 
mentally accessible. This critical particle radius could be 
further reduced either by increasing the fluid viscosity, 
or by using reactants molecules with a larger molecular 
weight [ll|. 

In order to further investigate the possibility for long- 
time self-propulsion, we need to go beyond the above 
linear stability analysis and check whether the asymp- 
totic state of the concentration field is compatible with 
swimming. In order to do so, numerical simulations of 
the full unsteady diffusiophoretic problem, Eqs. ©-(CO), 
are performed: the Stokes flow problem is solved explic- 
itly using the squirming mode decomposition, and the 
advection-diffusion problem is marched in time using a 
semi-explicit scheme, finite differences in the radial di- 
rections, and the Legendre spectral decomposition for 
the azimuthal dependence [311 ] . A small velocity per- 
turbation is imposed on the spherical particle initially at 
rest. Regardless of the amplitude of the initial perturba- 
tion, the system is stable and returns to its initial rest 
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FIG. 2: (Color online) Time-evolution of the instantaneous 
swimming velocity U(t) (thick solid line), and of three higher 
squirming modes, c n (r — l,t): mode n = 2, which corre- 
sponds to the magnitude of the induced stresslet (dashed line), 
n = 3 (dash-dotted line) and n — 4 (thin solid line). Results 
are displayed for Pe — 6 and AM = 1. The solute concen- 
tration is shown at three different times revealing the estab- 
lishment of a front-back asymmetry for an upward swimming 
motion. 



FIG. 3: (Color online) Long-time spontaneous swimming ve- 
locity, U°° (solid line), and magnitude of the long-time in- 
duced stresslet, E°° = — 4-7rMc2(r = l,t = oo) (dashed line), 
as a function of the Peclet number (with AM = 1). The 
steady state solute distribution around the particle is shown 
for Pe = 1, Pe = 5, and Pe = 15 for an upward swimming 
motion. The swimmer is always a "pusher" in the far field. 



state after perturbation for Pe < 4. In contrast, when 
Pe > 4, the swimming velocity of the particle grows ex- 
ponentially, and the growth rates obtained numerically 
are in quantitative agreement with our analytic predic- 
tions, see Fig. [TJ 

Turning to the long-time behavior, our computations 
confirm that the asymptotic state of the concentration 
field is compatible with locomotion. This is illustrated 
in Fig. [2] for Pe = 6, where we plot the time-evolution of 
the instantaneous particle swimming speed (thick solid 
line). A movie (supplementary information) and three 
typical snapshots in Fig. [2] illustrate the corresponding 
evolution of the solute concentration field. A steady 
state is reached and the swimming speed plateaus to a 
finite value in the long-time limit. Repeating the simu- 
lations for a range of Peclet numbers, we plot in Fig. [3] 
the nonlinear variation of the long-time swimming veloc- 
ity, U°° , as a function of Pe, demonstrating in particular 
the supercritical nature of the autophoretic instability. 
Note the non-monotonic variation with an optimal value 
of Pe sa 9 leading to the highest asymptotic swimming 
speed. For a given set of material parameters, and unlike 
autophoretic Janus particles [221 ] . the swimming speed 
here explicitly depends on the particle radius a on Pe. 
Physically, this difference arises because the concentra- 
tion gradient around an isotropic particle is not merely 
set by the particle size but instead by a dynamical scale. 

So far we focused exclusively on the impact of the con- 
centration field on the particle swimming speed. In order 
to address their collective dynamics, one would need to 



gain insight on the associated flow field away from the 
particle. The higher-order squirming modes are defined 
as c n (r,t) = (n+1/2) j_ 1 c(r, /i, t)L n (fx)dfx, where L n (/x) 
is the Legendre polynomial of order n. As discussed 
above, c\ is the (only) swimming mode. The flow field as- 
sociated with the mode C2 is the one with the slowest spa- 
tial decay, that of a stresslet (i.e. of a force dipole) [32j. 
Higher values of n correspond to higher-order flow sin- 
gularities in the far field. Repeating, for n > 2, the 
theoretical approach presented above for n = 1, we can 
obtain the unstable growth rates of each mode as a func- 
tion of Pe. The results reveal a hierarchy of supercritical 
instabilities corresponding to an infinite set of quantized 
critical Peclet numbers, Pe„ = 4(n+l), for mode n. This 
hierarchy of instabilities is illustrated in Fig. [T] where we 
plot the dependence of the growth rates for modes 2, 3, 
and 4 on the Peclet number. Notably, the stresslet mode 
becomes unstable at Pc2 = 12. 

Beyond linear analysis, the saturation of the swimming 
velocity results from the nonlinear evolution of concentra- 
tion fluctuations having a n — 1 symmetry into higher- 
order n-modes, as shown in Fig. [2j Although the Peclet 
number is below the critical value for all modes but the 
first one to be unstable (Pe = 6 in Fig. the nonlin- 
ear dynamics leads non-zero values for the other modes 
(modes 2, 3, and 4 are shown in Fig. [2]). In particular, the 
long-time value of the induced stresslet, is shown in 
Fig. [3] Independently of the signs of both A and M. , the 
unstable particle induces in the far field a flow with the 
symmetry of a "pusher" swimmer, similarly to flagellated 
bacteria [lj. 

Similarly to Marangoni flows which have been observed 
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for over 200 years to trigger the self-propulsion of cam- 
phor boats floating on water, in this letter we demon- 
strated that self-phoretic flows past isotropic particles 
show an instability to a spontaneous swimming state. 
The phenomenon discovered here could be exploited to 
design self-propelled "submarines" out of isotropic col- 
loidal particles. The instability mechanism is expected to 
be applicable generically to particles of different shapes 
or interacting with many chemical species, and suggests a 
simple experimental model system to carry out physical 
studies of active systems. 
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